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Single ion scattering contributions to the thermal conductivity of LiHoF 4 and 

LiY 96 %Ho4%F 4 

James NikkcQ and Brett EUmaotfl 

Kent State University, 105 Smith Laboratory, Kent, OH 
(Dated: February 2, 2008) 

We have performed extensive zero-field thermal conductivity measurements on single crystal sam- 
ples of LiHoF4 and LiY ge %Ho 4 %F4 below 2.3 K. By comparing these data to a single ion scat- 
tering model, we have shown that the thermal conductivity of LiYg 6 %Ho 4 %F4 is dominated by 
simple single-ion scattering while that of LiHoF4 shows additional contributions, possibly associ- 
ated with collective spin excitations. No entirely satisfactory model, however, is available to explain 
the thermal conductivity of the ferromagnet. 
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I. INTRODUCTION 

LiHoF 4 has been the focus of numerous experimen- 
tal and theoretical studies i£. The crystal field con- 
strains the angular momentum of the Ho 3+ ions to lie 
strictly along the c-axis. Furthermore, the large Ho-Ho 
spacing greatly reduces the exchange term, making the 
system an excellent approximation to a model dipolar 
Ising ferromagnet. The non-magnetic ion Y 3+ substitutes 
isostructurally and isoelectronically for Ho 3+ , leading to 
dilute 3-d Ising systems. The latter are, for low mag- 
netic ion concentrations, known to form a puzzling glassy 
system at low temperatures, whose signature is a nar- 
rowing distribution of barriers to relaxation as temper- 
ature is lowered. One goal of measuring thermal con- 
ductivity in these systems is to obtain coarse phonon 
spectroscopic data to probe the Ho 3+ ion energy lev- 
els in the low-temperature glass and ferromagnet. In- 
deed, the LiHoF 4 and LiY 96 %Ho4%F4 characteristic hy- 
pcrfine and magnetic energy level spectra are measured 
in the 100's of milli-Kclvins, making thermal phonons 
at temperatures ~ lif an experimentally simple way 
to access the characteristic energies of the system (sim- 
ple, that is, when compared to, e.g., millimeter pho- 
ton spectroscopy). The central questions we ask are (a) 
whether low-temperature collective excitations exist in 
the Ising ferromagnet LiHoF 4 below T c — 1.53-ftf and (b) 
whether the anomalous glassy behavior of the dilute sys- 
tem LiY g6 %Ho4%F4is due to rearrangements of correlated 
spin clusters. In order to address these issues, we must 
analyze the data in such a way as to deconvolve the rather 
complex single Ho 3+ ion contribution to the thermal re- 
sistance. We will accomplish this via a detailed analysis 
of the crystal field splittings and resulting single-ion ex- 
citations that cause phonon scattering. 



k = J C(k)u(k)Z(k) d 3 k. 



(1) 



Our initial approximations are to assume that the system 
is structurally isotropic with respect to k and that only 
acoustic phonons with a linear dispersion relation con- 
tribute at the experimental temperatures. We can then 
calculate the thermal conductivity in terms of an integral 
over the phonon frequency, u>, as 
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The specific heat, C{lj) contribution from phonons 
with energy between cu and u> + du, is the derivative of 
the mean thermal energy density of the lattice 3 , 



C(ui) du> 
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Using equation [2 the thermal conductivity is now 
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Our task is to find l(u>), taking into account both 
single-ion magnetic scattering and scattering due to, e.g., 
crystalline defects. For reference, we reproduce in tablelTTl 
a synopsis 3 of the u> dependence of various structural 
scattering mechanisms. 

We also assume the correctness of Mathiesen's rule, 
i.e., that contributions the mean free path are added re- 
ciprocally. 



II. SCATTERING MODEL 

We start with the general form of the thermal conduc- 
tivity, k, in terms of the phonon heat capacity, C(k), the 
mean free path, Z(k), and the the speed of sound, u(k), 
written as an integral over the phonon wave vector, k: 



III. DETAILED SINGLE ION CALCULATION 

The thermal resistance of an isolated Ho 3+ ion may 
be numerically estimated from the electronic transition 
probabilities of the ions when they are excited by thermal 
phonons. We will not attempt here to calculate from first 
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Type of scattering 


l(u>) oc 


k(T) <x 


External boundaries 


const. 


rj-i3 


Grain boundaries 


const. 


rj-i3 


Stacking faults 


oj- 2 


T 


Conduction electrons 


Id 




Point defects 


0J~ 4 


1/T 


Umklapp processes 


OJ 


T 3 e a/T 



TABLE I: Table of the functional form of the mean free path, 
I, on oj. 



principles the electron-phonon coupling constants. In- 
stead, we will use a semi-phenomenological model to as- 
certain whether a single-ion scattering term can account 
for the observed thermal conductivity. The 5 2g multiplet 
of the Ho 3+ ions in LiHoF4 is comprised of a ground 
state doublet followed by two singlets at about 9 and 
39 K. The energy levels above these are of sufficiently 
high energy that we can neglect them at the experimental 
temperatures. The ground state doublets are each split 
into 8 hyperfine sublevels which in turn may be split rel- 
ative to each other due to a Zeeman interaction with an 
external or internal magnetic field. The magnitude of the 
ground state magnetic splitting in the ferromagnet was 
measured by Battison et al. 4 . Using an optical spectro- 
scopic technique, they found it to be proportional to the 
magnetization (i.e., mean field) with a maximum at zero 
temperature equal to 2.6 cm~ l , or 3.7 K. 

The energy levels are calculated by solving the eigen- 
value problem for the 5 /g single ion Hamiltonian. Wc 
follow the method described by Giraud et al&. In the 
\J,m, I,n) basis, the Hamiltonian can be approximated 
as the sum of the crystal field, Zeeman, and hyperfine 
contributions, 



H — HcF + Hzee + Hhf- 



(•5) 



The crystal field contribution is written in term of the 
Stevens' operators, 0™, and the crystal field parameters, 
B™ ■ The resulting crystal field Hamiltonian is 
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(6) 



where a, (3, and 7 are ion specific coefficients calculated 
by Stevens 6 . For Ho 3+ , a = -1/450, [3 = -1/30030, 
and 7 = —1/3864861. (The general form for a and (3 
can be found on page 253 of Hutchingsi.) The crystal 
field parameters, B™ , are given by Gifeisman et al. in 
reference^ and were measured through high resolution 
optical spectroscopy. Their values are: B\ = 273.9 K, 
B% = -97.7 K, Bl = -6.5 K, B\ = -1289.1 K, and 
B\ = -631.6 K. 
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in I J, m) basis 



3J Z 2 - J(J + 1) 
35J Z 4 - 30J( J + 1) Jl + 25J Z 2 
-6J(J+1)+3J 2 (J + 1) 2 

231J Z 6 -315J(J + 1)J 4 

+ 735J 4 + 105J 2 (J + l) 2 Jl 

- 525J(J + 1)J 2 +294J 2 

-5J 3 (J+1) 3 + 40J 2 (J+1) 2 

-60J(J + 1) 
\ [(11 Jl - J(.J + 1) - 38) (j+ + Jl) 

+ (4 + Ji) (11J 2 - J(J + 1) - 38)1 



TABLE II: Table of the Stevens' operators and equivalents. 



The Stevens' operators, OJ" are given in reference |y 
and and the ones of interest here are summarized in 
tabic ITTT1 

The Zeeman contribution to the Hamiltonian is simply 



Hze 



(7) 



where the g-factor gj = 5/4 for the J = 8 ground state, 
fiB is the Bohr magneton, and H is the magnetic field 
due to applied external fields and the internal field in the 
ferromagnetic state. Due to the symmetry of the crystal, 
the internal magnetic field due to ferromagnetic order is 
directed along the (Ising) c — axis. Therefore, in the 
absence of an external field, and within the mean-field 
approximation, we may write 



H Z ee = -g,jiLB^-rM{T)J z 

M 



(8) 



where Mq is the saturation magnetization at zero tem- 
perature, Hq is the effective maximum internal field at 
zero temperature, and M(T) is the temperature depen- 
dent magnetization. Note that this is the only term in 
the Hamiltonian where the temperature enters into the 
calculation of the energy spectrum. (The thermal broad- 
ening of energy levels described below also depends on T.) 
The magnetization of LiHoF4 is calculated from a fit to 
the magnetization data measure in reference U The zero 
temperature internal field, Hq, is approximately 0.33 T 
which gives the correct splitting of 3.7 K at zero temper- 
ature as measured by Battison et al. 0. Note that since 
our sample geometry is a long bar along the c-axis, we 
do not apply a demagnetization correction. 
The hyperfine contribution may be written 



H, 



AjJ-I. 



(9) 



The coupling parameter Aj was found to be approxi- 
mately 39 mK by Mennenga et al*& for LiHoF4 from 
specific heat measurements. Magarino et ali^ found ap- 
proximately the same value for 2% Ho 3+ diluted into 
LiYF4 using electron paramagnetic resonance. 
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For the Ho 3+ ions in LiHoF4 and LiY g6 %Ho4%F4, 
J = 8 and / = 7/2. Therefore, the Hamiltonian takes 
the form of a 136 x 136 Hermitian matrix. Numerical di- 
agonalization gives the respective eigenvectors, which in 
turn allow the calculation of pertinent expectation val- 
ues. Specifically, it is important to know the expectation 
value of the z-component of the nuclear angular momen- 
tum, (I z ), because transitions that change (I z ) a large 
amount are strongly suppressed. Figure ^ is a plot of the 
ground state doublet and the first excited state eigen- 
values plotted against the expectation values (J z ). The 
value of (I z ) is written next to each point. The values for 
(J z ) are uniformly half integer as we would expect. Note 
that the scale of the (J z ) axis on the top of the plot is one 
tenth that of the bottom to show the form of the excited 
state more clearly. Figure |3 is similar to figure ^ but 
with a magnetic field turned on in the z direction. The 
field used here is about 0.3 T, equal to that which would 
be generate internally by LiHoF4 at zero temperature. 
As we will discuss, only transitions between energy levels 
with A J = will be included in our computations. 

There is one important note here about the energy 
eigenvalues in our particular case. The crystal field pa- 
rameters measured by Gifeisman et al& were measured 
through high resolution optical spectroscopy and they 
did not observe the 9.5 K excited state in their lumines- 
cence spectra. Their parameters inserted into the calcu- 
lation therfore predict that this unconstrained level sits 
at 14.5 K above the ground state doublet. Since this 
is the only state so effected, we adjust our eigenvalues 
by subtracting 5 K from this state before using them in 
further calculations. 
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FIG. 1: Plot of the energy eigenvalues arranged in terms of 
the (J z ) expectation values. The labels on the points are the 
(I z ) expectation values. The bottom plot shows the ground 
state doublet, and the top plot is the I s * excited state. 
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FIG. 2: Similar to figureQbut with magnetic splitting equal 
to what would occur in LiHoF4 at K. Note that the scale 
of the ( J z ) axis is the same in the top and bottom of the plot. 



assume that there are only three scattering processes, 
giving us three independent Us to add. The first is due 
to phonon scattering from the sample boundaries. This 
term is simply the geometrically averaged sample size, 
Ibr, and is independent of w. The second term is from 
scattering off of point defects, l p d, which gives a term 
proportional to 1/cu 4 (see table [nj. As we shall see, this 
term is necessary to obtain agreement with the high tem- 
perature range of our data (though higher dimensional 
defects might also suffice, we assume point defects for 
simplicity). The last process is due to the non-trivial 
scattering from the Ho 3+ ions, lion- Therefore, 



I = 



1 



l 

lpd 



(10) 



We choose to write the mean free path due to the single 
ion scattering process as li on = l m in/ Pscatt- The scale 
factor l m i n is similar to the minimum distance that a 
phonon can go without running into a Ho 3+ ion, and 
Pscatt is the probability that the phonon will scatter off 
of that ion. We define l p d = -^r, and will use £ as a 
fitting parameter. Our complete expression for the mean 
free path is then 



P scatt (uj,T) 



(11) 



After computing the energy level spectrum at a par- 
ticular temperature, we calculate the partition function 
which is obtained by simply summing over each level, Ei, 
integrated over the level density, Di(e), 



Recall that the goal is to calculate a temperature de- 
pendent mean free path, 1(lo), that is to be inserted in 
to the integral in equation 0] In the present case, we 



Z 



E 



(12) 
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We assume that the individual energy levels are ther- 
mally broadened (due to crystal field fluctuations from 
thermal motion of the atomic neighbors of the Ho 3+ ions 
in the lattice) with a Gaussian energy distribution. For 
simplicity wc further assume that all of the levels are 
have the same width, AE, resulting in the following en- 
ergy density function, 



A(e) 



1 



ae^h 

Our partition function now becomes 



(13) 



(14) 



The scattering probability, P SC att(Ep,T), is a function 
of the phonon energy, Ep — Hlo, and the temperature of 
the system, T. We calculate it by performing a double 
sum over all the energy levels, calculating the scattering 
probabilities between each pair. 

P scatt (E P ,T)=s(\{lW)-{lW)\) (15) 

r p oo 

X E / e-*I>y(e)(fe 
i,j lJ -°° 

/OO 
-oo 

where D i3 (f) = Di(e) Dj(e + Ep) and Dji(e) — 
Dj(e) Di(e — E p ). This particular form takes into ac- 
count stimulated emission from the higher level to the 
lower one as well as the direct excitation of the ion. The 
delta function allows for excitations only between levels 
that have the same nuclear angular momentum expecta- 
tion values (see Fig's ^ and[2J- Slightly decreasing this 



selectivity, say by allowing 



to vary by one 



or two, makes little difference to the current model ex- 
cept at the lowest temperatures (< 100 mK). Using the 
level density mentioned above, the scattering probability 
is found to be 



P scatt (E P ,T) = ^ 2AE2 
x < e^ 2 - e-^/ 2 



(16) 



| e /3Bp/2 _ e -/3B P /2| 
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The energy level width, AE, requires attention here. 
The thermal energy due to deformation of the lattice is 
essentially 



E, 



k h T 



-C e 2 



(17) 



where Co is an averaged elastic constant, and e is an 
averaged strain. For small lattice deformations changes 
in the crystal field energy are therefore proportional to 
the strain, 



SEcf oc e. 



(18) 



We collect the constants, and designate the thermal 
broadening as 



AE = AVT. 
with A as a fitting parameter. 



(19) 
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FIG. 3: Plots of the contributions to the heat capacity calcu- 
lation for T = 2 K. The various panels are described in the 
text. 

The thermal conductivity is calculated at each tem- 
perature by the following procedure. First the vari- 
ous temperature dependent parameters (internal field for 
LiHoF4and the broadening parameter A) are set and the 
Hamiltonian eigenvalue problem is solved. The integral 
in equation 0] is computed by inserting the mean free 
path, l{u>), from equation 1111 The scattering probabil- 
ity is then calculated at each point in the integral using 
ea uat ion 1161 where Ep = fru>. 

We have used an integration cut-off of u> — lQksT/fi, 
which is quite sufficient as seen in Figure |21 Panel a is 
the heat capacity part of the integrand which essentially 
vanishes for higher frequencies. Panel b is the scattering 
probability from equation 1161 and panel c is 1(lo) from 
equation II II The last panel, d, is the integrand to equa- 
tion 01 so the area of the filled portion represents the 
total thermal conductivity at this temperature. The ef- 
fect of transitions between Ho 3+ levels is quite apparent 
in the mean free path data, effectively removing specific 
frequency bands from the thermal conductivity integral. 
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IV. EXPERIMENTAL DETAILS 

Thermal conductivity was measured by the standard 
gradient method using a heater and two thermometers. 
Our heaters were 200 ft metal film on Kapton strain 
gauges and the thermometers used were 1 kVl to 10 kfl 
R11O2 surface mount resistors, depending on the tem- 
perature range of interest. Copper bars approximately 
150 fim thick were silver epoxied to the sample divid- 
ing it approximately into thirds. These bars defined the 
distance over which the temperature gradient, AT = 
T2 — Ti, was measured. The heaters were wired with 
superconducting NbTi wires to eliminate ohmic heat- 
ing and to minimize heat leaks, while thin Evanohm 
wires were used on the thermometers. To check for 
thermal leakage through the wires, the measurement of 
LiY g6 o/ Ho 4 %F4 was repeated with the wires doubled 
in length. The results of the thermal conductivity be- 
tween the two sets of measurements were within the noise 
of either experiment. Our two samples (pure and di- 
lute) were optically oriented and cut such that the long 
length was along the c — axis. The pure LiHoF 4 was 
1.46 x 1.56 x 5.4 mm 3 . This gives a cross section of 
0.0228 cm 2 . The distance, L, between the thermometers 
was 0.22 cm. The 11.7 mg dilute LiYg 6 %Ho 4 %F4 sample 
was 0.93 x 0.75 x 4.7 mm 3 , but slightly trapezoidal. The 
effective cross section was calculated by using measure- 
ments from a digital photograph of the sample. The final 
value used was A = 0.0066 cm 2 . The distance between 
thermometers for this sample was also 0.22 cm. Most ex- 
periments were conducted in a 3 He cryostat. Data below 
?>10mK were obtained in a helium dilution refrigerator. 



V. EXPERIMENTAL RESULTS 



A. Pure LiHoF 4 




1 1.5 2 

Temperature (K) 

FIG. 4: Measured thermal conductivity, k, for LiHoF4. The 
bottom panel is a magnification of the top to better illustrate 
the structure below the Curie Temperature, Tc- 



The accumulated data for a number of experiments 
are plotted in figure The top panel is the entire range 
of data obtained, and the bottom is a magnification of 
data below 2 K. One can easily see the ferromagnetic 
transition on this plot at ps 1.54 K. One interpretation 
of the overall structure of the data is as follows; as we go 
down in temperature, the thermal conductivity decreases 
due to the T 3 term in the specific heat. At this point, the 
vast majority of phonon-driven transitions are between 
the hyperfine ground state sublevels. When we reach 
Tc, the ground states start Zeeman splitting, and there 
are a number of new transition possibilities opened up, 
so the conductivity drops rapidly. Note that the peak in 
C(uj) in figure |31 occurs at approximately u> = AkbT/h. 
This implies that at Tc, about 1.5 K, the majority of 
phonons have energies around 6 K. However, just above 
Tc, the Ho 3+ energy spectrum looks like that shown in 
figure m where there are no scattering possibilities for 
phonons between about 2 K and about 8.5 K. Below Tc, 
the excluded temperature region quickly shrinks due to 
the splitting, and simultaneously, the peak of the phonon 
spectrum decreases. At about 1.2 K, the ground states 
have split enough that transitions are suppressed and the 
conductivity starts increasing. Eventually, this rise is 
expected to be countered by the T 3 factor from having a 
sample of finite size. While this qualitative explanation 
was the primary motivation to develop a single ion model 
of the thermal conductivity we shall see that this picture 
is not sufficient to explain all of the data. 




1 2 
Temperature (K) 

FIG. 5: log- log plot of k for LiHoF4 along with some lines of 
T n for comparison with table HTl 



Figure [5] shows the whole range of data on a log-log 
plot, along with lines proportional to T, T 2 , and T 3 . Re- 
call from table [n] the various structural scattering pro- 
cesses that can contribute to the thermal conductivity. 
As one can see, the data do not follow any one of these 
power laws over any significant temperature range. We 
did not take measurements of the pure LiHoF4in the di- 
lution refrigerator, so we can not see the approach to T 3 
in the low temperature limit. 
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Temperature (K) 

FIG. 6: Measured k for LiY g6 %Ho4%F4. The bottom panel is 
a magnified view of the top. 



B. Dilute LiY g6 %Ho4%F4 

We have again plotted our data in two temperature 
ranges in figure El This plot is not very illuminating 
by itself since it goes so quickly and smoothly to zero. 
There is no ferromagnetic transition, so there are no 
sharp structures as there are for LiHoF 4 . (The gap in 
the data between 1.3 K and 1.5 K is due to T-control is- 
sues in the 3 He cryostat.) As in the previous section, the 
data becomes increasingly noisy at higher temperatures 
due to the reduced sensitivity of the thermometers and 
the rapidly increasing thermal conductivity. Note that at 
2 K the thermal conductivity of LiYgg%Ho 4 %F4 is about 
10 times that of L1H0F4, but slightly smaller at 0.5 K. 
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FIG. 7: log-log plot of k for LiYgg%Ho 4 %F4 along with some 
T n lines for comparison. 

Figured is a plot of log(«;) versus log(T) to check for 
regions where the thermal conductivity is dominated by 
a single power law. It appears that at the lowest tem- 
peratures the data may be approaching T 3 , but going to 
lower temperatures would be necessary to confirm this. 



More importantly, the log plot brings out structures that 
must be reproduced by our single-ion model if it is to 
describe the data. 

VI. ANALYSIS 

Our goal is to find a set of parameters that provide the 
best fit for both our LiHoF4and LiY 96 %Ho 4 %F4 data si- 
multaneously. Fitting parameters include £, which deter- 
mines the contribution from point defect scattering, and 
Imin, which determines the size of the contribution from 
the Ho 3+ ion scattering. P sca tt also contains a parameter 
that can be varied, as does the the thermal broadening, 
Eq. Since A is a single ion property, it should be 

independent of the concentration of the Ho 3+ ions. Sim- 
ilarly, Imin should be proportional to the inverse of the 
cube root of the concentration, fixing its relative value. 
The effective mean free path due to the finite size of the 
sample, Z& r , is calculated from the sample dimensions and 
is fixed. We take Aj from equation [5] to be a constant, 
38.6 mK , since this value has been measured and pub- 
lished using at least two different methodsiSiii. As men- 
tioned, Hq 0.33T from equation [S] was indirectly, and 
somewhat imprecisely, measured by Battison et alX. 

We will first discuss the dilute glassy system 
LiY 96% Ho 4 %F 4 . 

A. Dilute LiY g6% Ho 4 %F4 




Temperature (K) 



FIG. 8: Plot of nv 2 /aT 3 and best fit for LiY 96% Ho 4 %F 4 . Inset 
is K. 

Figure [SJ is a plot of k and KV 2 /aT 3 for 
LiY 96 %Ho4%F 4 along with the best model fit. Here, 



which gives us k/T 3 roughly in units of the mean free 
path (the equivalence would be exact of I did not de- 
pend on lo). As we can see, the fit is quite good with 
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this parameter set. The best fit parameters are A — 
0A5±0.03K 1 / 2 , l mm = 38±1 fm, and £_ = 15±2 m-K 4 . 
The uncertainties were estimated by varying the individ- 
ual parameters to obtain the same minimum % 2 , and do 
not reflect any correlations between the parameters. As 
we can see, the scattering due to point defects is a fairly 
small contribution. For example, a 1 if phonon will have 
a mean free path, due to this process, of about 15 me- 
ters. Only for phonons with energies greater than about 
12 K does the mean free path to become smaller than the 
physical size of the sample. Because of this, the initial 
fitting was done with £ = oo, after which £ was adjusted 
to minimize the x 2 . 




Temperature (K) 



FIG. 9: Plot of KV 2 /aT 3 and variations on the best fit for 
LiY 96 %Ho 4 %F4. The parameters are described in the text. 

Figure demonstrates how each parameter effects the 
model. Unless noted, all lines shown here use the "best 
fit" parameters. The dashed line has the point defect 
scattering turned off by setting £ = oo. Note that this 
line follows the best fit line fairly well for temperatures 
below the peak at 900 mK, and is identical for temper- 
atures less than 500 mK. The dot-dashed line demon- 
strates the effect of turning off the Ho 3+ single ion scat- 
tering process by by setting l m i n — oo. If we turned off 
the point defect scattering at the same time, we would 
obtain a horizontal line at kv 2 / aT 3 =1(^=1 A9 mm. From 
this, we see that the single ion process alone is responsible 
for the structure of the dip at 270 mK and the addition 
of the point defect scattering causes a turn-over in the 
data to produce the peak at 900 mK. 

In the very low temperature region, less than 80 mK, 
the model demonstrates a flattening behavior that is not 
present in the data. This is due to the value of lb r , since 
it sets the overall maximum mean free path. Interest- 
ingly, increasing lf, r beyond its physical value does not 
improve the overall model, as this effectively scales the 
entire fit. There currently is no parameter set in with 
this model that reproduces the T < 80 mK region while 
retaining agreement with the rest of the data set. One 
explanation for this could be related to the surface con- 
dition of the sample. At very low temperatures, phonons 



will have an increasing wavelength, which will eventually 
become larger than the surface roughness. When this 
occurs, the boundary scattering becomes specular and 
the effective mean free path due to the finite size of the 
sample increases. This would result in a rapid upturn 
of k/T 3 when this limit is reached, consistent with our 
measurements. 

We believe that this single ion model describes the 
thermal conductivity of LiYg 6 %Ho4%F4 fairly well be- 
low 2.5 K, and demonstrates the relative importance of 
scattering from point defects. We see no sign of the "clus- 
ter glass" proposed by Reich et. af •. They observed 
narrowing in the magnetic susceptibility with decreasing 
temperatures between 300 mK and 150 mK, but we see 
no signs of collective behavior in this temperature range. 



B. Pure LiHoF 4 
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FIG. 10: Plot of Kv 2 /aT 3 and three fits for LiHoF 4 . Param- 
eters are described in the text. 

We now move to the ferromagnet LiHoF4. Because 
A should be independent of the concentration of the 
Ho 3+ ions, and l m in should be proportional to the in- 
verse of the cube root of the concentration, we will start 
with A = 0.45 K 1 / 2 , l mm = 3 VWo ■ 38 = 13 fim, and 
£ = 15 m ■ K 4 . from the analysis of LiY 96 %Ho4%F4. We 
also turn on the Zeeman splitting by setting H from 
equation to 0.33 T. Figure [TTJ is a plot of nv 2 /aT 3 
versus temperature for LiHoF4 with three lines from 
our model. Note that the vertical axis is on a log scale 
so that the structure of the ferromagnetic transition can 
be seen. The solid line is our "best fit" with the above 
parameters. The dashed line uses the same parameters 
except for £ = oo to see the effect of turning off the point 
defect scattering. The dot-dashed line has l m in = oo 
to demonstrate the effect of turning off the Ho 3+ sin- 
gle ion scattering contribution (the individual contribu- 
tions are more profitably viewed using the color figure 
available online). As we can clearly see, the "best fit" 
is poor. Furthermore, no amount of adjustment of the 
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parameters appears to yield significantly better agree- 
ment. We do, however, see some qualitative successes 
with this model of L1H0F4. The behavior immediately 
below the Curie temperature is generally what we pre- 
dicted from section lV Al and we see a reasonably accurate 
roll off in k/T 3 at the lowest temperatures. One possi- 
bility is that there are one or more additional scattering 
mechanisms needed in the model. For example, below 
the Curie temperature, domain formation could intro- 
duce a temperature dependent grain boundary-like scat- 
tering process. However, this would have no effect above 
Tc, where the deviation from the fit is most pronounced. 
(Of course, it is possible that difficulties present below 
Tc in the model result in best-fit deviations above the 
Curie point.) The second, and potentially most interest- 
ing possibility is that phonon-excited spatially correlated 
spin-flips account for the enhanced scattering. Note that, 
unlike isotropic or nearly isotropic systems, conventional 



spin-waves are not expected in an ensemble of Ising spins. 
Indeed, the nature of (potentially collective) magnetic 
excitations above and below Tc in this dipolar-coupled 
system are unclear, as is the question of whether such 
magnetic scattering could account for the behavior we 
observe. Again, one would naively believe such effects 
to strongly dominate below Tc- At higher temperature 
it is unclear what multi-spin excitations are available to 
couple to the phonon current. 

In conclusion, a detailed semi-phenomenological model 
of single-ion phonon scattering processes is in excel- 
lent agreement with thermal conductivity data on di- 
lute LiYg 6 o/ Ho 4 °/ F4, with no evident signature of low- 
energy collective excitations involving random spin clus- 
ters. Similar data on the pure ferromagnetic system 
LiHoF4, however, are not in quantitative agreement with 
the independent ion theory and will presumably require 
a study of collective effects in this dipolar Ising system. 
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